Int. J. Mol. Sci. 2014, 15, 9242-9254; doi:10.3390/ijmsl5069242 



OPEN ACCESS 



International Journal of 

Molecular Sciences 

ISSN 1422-0067 

www.mdpi.com/journal/ijms 

Article 

Genetic Isolation among the Northwestern, Southwestern and 
Central-Eastern Indian Ocean Populations of the Pronghorn 
Spiny Lobster Panulirus penicillatus 

Muhamad Fadry Abdullah *'*, Alimuddin 2 , Mohamed Muththalib 3 + , Adnan Jameel Salama 4 
and Hideyuki Imai 5 

1 Graduate School of Engineering and Science, University of the Ryukyus, Nishihara, 
Okinawa 903-0213, Japan 

Faculty of Fisheries and Marine Science, Bogor Agricultural University, Bogor 16680, Indonesia; 

E-Mail: alimuddin_alsani@yahoo . com 

Ministry of Fisheries and Agriculture, Male 20125, Maldives; 

E-Mail: mohamedmuththalib@gmail.com 

4 Faculty of Marine Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia; 
E-Mail: ajsalama@kau.edu.sa 

5 Faculty of Science, University of the Ryukyus, Nishihara, Okinawa 903-0213, Japan; 
E-Mail: imai@sci.u-ryukyu.ac.jp 

t Present address: Kooddoo Fisheries Maldives Ltd., Block 389 Ujaalaa Higun, Male, Maldives. 

* Author to whom correspondence should be addressed; E-Mail: fadrynankeren@yahoo.com; 
Tel.: +81-98-895-8554; Fax: +81-98-895-8576. 

Received: 10 April 2014; in revised form: 14 May 2014 /Accepted: 16 May 2014 / 
Published: 26 May 2014 

Abstract: The pronghorn spiny lobster Panulirus penicillatus is a highly valuable 
species which is widely distributed in Indo-West Pacific and Eastern Pacific regions. 
Mitochondrial DNA control region sequences (566-571 bp) were determined to investigate 
the population genetic structure of this species in the Indian Ocean. In total, 236 adult 
individuals of Panulirus penicillatus were collected from five locations in the Indian Ocean 
region. Almost all individuals had a unique hap lo type. Intrapopulation haplotype (h) and 
nucleotide (n) diversities were high for each locality, ranging from h = 0.9986-1.0000 and 
7T = 0.031593-0.043441. We observed distinct genetic isolation of population located at the 
northwestern and southwestern edge of the species range. Gene flow was found within 
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localities in the central and eastern region of the Indian Ocean, probably resulting from an 
extended planktonic larval stage and prevailing ocean currents. 

Keywords: pronghorn spiny lobster; Panulirus penicillatus; genetic diversity; population 
structure; Indian Ocean 



1. Introduction 

The pronghorn spiny lobster Panulirus penicillatus is widely distributed in the Indo-West Pacific 
and Eastern Pacific regions. It is found in tropical and adjacent regions from South-Eastern Africa, 
the Red Sea, Southern India, the Southeast Asian Archipelago, Japan, Northern Australia, and the 
Southern and Western Pacific Islands, to Hawaii, the Galapagos Islands, and other islands of the Eastern 
Pacific [1-3]. P. penicillatus supports considerable fisheries in Indo-Pacific regions. The species is 
fished wherever it occurs. Much attention has been paid to biological investigation of young and adult 
stages. Recognition of stocks is complicated by the life cycle of pronghorn spiny lobster, primarily 
because of the potential for long-range dispersal during its planktonic larvae. The early life history of 
spiny lobsters consists of a drifting larval period adapted for a relatively long-term stay in the open 
ocean, extending from several months to more than a year, with many possibilities for dispersal 
through ocean currents [4,5]. 

Several molecular techniques, including allozymes, restriction fragment length polymorphism 
(RFLP), mitochondrial DNA sequencing, and microsatellites, have been used in studying genetic 
variability in marine organisms at the population level [6,7]. Mitochondrial DNA sequencing, 
particularly of the most rapidly evolving and highly variable control region, has been a useful tool 
for population genetic studies of many terrestrial and aquatic organisms [8]. The highly polymorphic 
mitochondrial control region sequence has been used previously for population genetics studies in 
spiny lobster species [9-12]. 

There recently been studies of P. penicillatus, East- West Pacific differentiation using 16s and COI 
(cytochrome oxidase I) marker was described by Chow et al. [13], and this study was corroborated 
with the study of genetic population structure of P. penicillatus in the Pacific region using highly 
polymorphic marker, mitochondrial control region by Abdullah et al. [3]. In this study, we focused on 
the Indian Ocean region off. penicillatus. 

There have been no previous population genetic studies of P. penicillatus in the Indian Ocean 
region. The aim of the present study is to investigate population genetic structuring in a widely 
distributed P. penicillatus in the Indian Ocean region. We address this question using mitochondrial 
control region sequences data from Red Sea, Jeddah, Saudi Arabia; Fort Dauphin, Madagascar; 
Nilandhoo atoll, Maldives; Aceh and Java Sea, Indonesia (Figure 1). Given the economic importance 
of P. penicillatus, this genetic information may be valuable for long-term fisheries' management decisions. 
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Figure 1. Map of the Indian Ocean (left): showing SEC — South Equatorial Current, 
SECC— South Equatorial Countercurrent, SEMC (NEMC)— Southeast (Northeast) 
Madagascar Current, NMC (SMC) — Northeast (Southwest) Monsoon Current, JC — South 
Java Current, ocean current was drawn referring to Schott and McCreary [14]. 
Panulirus penicillatus haplotype network (right); Sizes of the circles are proportional 
to the frequency of each haplotype; lengths of the lines are relative to the number of 
mutations between haplotypes. 
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2. Results and Discussion 



2.1. Genetic Diversity Analysis 

A total of 566-571 bp sequence data were obtained from 236 specimens of P. penicillatus. 
Of the 577 aligned base pairs, 246 variable sites (including indels) were found: 172 were parsimony 
informative and 74 were singletons. In total, 231 haplotypes were identified, and nearly all sequenced 
individuals had a unique haplotype (Table 1). The nucleotide sequences of the haplotypes were 
deposited in the DNA Data Bank of Japan (DDBJ) (accession numbers AB922639-AB922816, 
AB689451-AB689503). 

As expected, this fragment is AT -rich in nucleotide composition with ratios of 37.8% adenine, 
31.8% thymine, 19.4% cytosine, and 11% guanine, consistent with previous reports of an AT-rich 
control region of the mitochondrial genome in many invertebrates, including crustaceans [7,15]. 

The haplotype network showed a geographic structure with two main groups (Figure 1): (1) the 
northeastern edge of the population; Red Sea cluster and (2) other localities from the southwestern to 
eastern Indian Ocean (Madagascar, Maldives, Aceh and Java) that are closely together. 
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Table 1. Sample abbreviation (ID), sample locality, sample size (N), year of collection, numbers of haplotypes, haplotype diversity (h), 
nucleotide diversity (n), Tajima's D and Fu's F Test for Panulirus penicillatus mitochondrial DNA control region. 



ID 


Locality 


N 


Year 


No. of Haplotypes 
(Unique) 


/i±SD 


7T±SD 


Tajima's D 
(p-Value) 


Fu'sF 
(p-Value) 


JD 


Red Sea, Jeddah, Saudi Arabia 


46 


2011 


45(45) 


0.9990 ± 0.0048 


0.031593 ±0.015851 


-0.79131 (0.215) 


-24.20576 (0.000) 


MG 


Fort Dauphin, Madagascar 


49 


2014 


49(49) 


1.0000 ±0.0041 


0.043441 ± 0.021533 


-0.97966 (0.168) 


-24.14011 (0.000) 


MD 


Nilandhoo atoll, Maldives 


38 


2013 


37(36) 


0.9986 ± 0.0065 


0.034957 ±0.017561 


-1.05781 (0.127) 


-19.70407 (0.000) 


AC 


Aceh, Indonesia 


48 


2011 


47(47) 


0.9991 ±0.0045 


0.038100 ±0.018974 


-1.17981 (0.098) 


-24.14392 (0.000) 


JV 


Java Sea, Indonesia 


55 


2008 


53(52) 


0.9993 ± 0.0036 


0.038855 ±0.019284 


-1.01996 (0.153) 


-24.12306 (0.000) 




Total 


236 
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Genetic diversity, as shown by haplotype diversity, was found to be very high in P. penicillatus. 
The values generated in this study are comparable with those derived from the mtDNA control region 
sequences for the crustaceans Farfantepenaeus aztecus, Litopenaeus setiferus and Sergia lucens [16,17]. 
Haplotype diversity (h) within the geographical populations was high, ranging from 0.9986 in 
Nilandhoo atoll, Maldives to 1.0000 in Fort Dauphin, Madagascar. Nucleotide diversity (n) was 
generally high, ranging from 0.031593 in Red Sea, Jeddah, Saudi Arabia to 0.043441 in Fort Dauphin, 
Madagascar (Table 1). There were few shared haplotypes among locations belonging to the 
central-eastern Indian Ocean. 

In the present study, the mtDNA control region exhibited high genetic variability (h = 0.9986-1.0000 
and 7r = 0.031593-0.043441), which are relatively high values compared with the control region 
diversities of other Panulirus spp. (P. inflatus: h = 0.957-1.0000 and n = 0.019-0.028 [10]; 
P. argus: h = 0.667-1.0000 and n = 0.003-0.066 [7], P. homarus: h = 0.985-1.0000 and 
71 = 0.031-0.041 [12]) and similar to Pacific region populations of P. penicillatus (h = 0.9996-1.0000 
and 7T = 0.0255-0.0448) [3]. The control region is a non-coding region that includes signals necessary 
for replication of a molecule and is the most rapidly evolving region of mtDNA [18]. 

The high genetic variability of the control region in P. penicillatus may be maintained by its large 
population size. In fact, because nearly all individuals represented a unique haplotype, an investigation 
based on sequencing the control region would require a very large sample size to detect common 
haplotypes and associate them with particular geographic areas [10]. The evolutionary pattern of 
interhaplotypic divergence will almost always produce a biogeographic signature, and an analysis 
that recognizes genetic distance among alleles will provide additional resolution, especially for 
hypervariable markers [19]. 

2.2. Population Structure Analysis 

Genetic analysis of mitochondrial control region revealed distinct population at the northwestern 
and southwestern margin of geographic range of P. penicillatus with other localities. This suggests that 
the geographic isolation may be due to contemporary isolation/gene flow or prior isolation with recent 
gene flow, although the overall level of genetic exchange must be below that required to homogenize 
the populations. 

Congruent with the pattern observed in the network haplotypes (Figure 1), we found strong regional 
structuring between Red Sea and other Indian Ocean populations of P. penicillatus, which explained 
high genetic variations (51%) between the regions. High genetic variation (85%) was also observed in 
previous study between eastern Pacific and west-central Pacific population of P. penicillatus [3]. 
To analyze the sub-structure in the other area we attempted three types of regional grouping: 
(1) Madagascar and Maldives/ Aceh/Java; (2) Madagascar/Maldives and Aceh/Java; and 
(3) Madagascar/ Aceh and Maldives/Java. Regarding these groupings, genetic variation was relatively 
greater between Madagascar and Maldives/ Aceh/Java (2.41%) than the variations between other 
combinations (%var -0.49 and -0.75) (Table 2). 
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Table 2. AMOVA results showing degree of freedom (df), variance components (Var.), 
percent variation (% Var.) and (^-statistics for Panulirus penicillatus (1000 permutations). 



Source of Variation 


df Var. % Var. 


^-Statistics 


P 


Value 


All sites (Red Sea/Madagascar-Maldives-Aceh-Java) 


Among groups 


1 11.50819 51.31 




0.51313 


P = 


0.19746 


Among populations 


3 0.18886 0.84 


®sc = 


0.01730 


P = 


0.00098 


Within populations 


231 10.73030 47.84 


@ST = 


0.51313 


P = 


0.00000 


Substructure grouping test: (Madagascar-Maldives/Aceh/Java) 


Among groups 


1 0.27591 2.41 


0 CT = 


0.02408 


P = 


0.24047 


Among populations 


2 0.03850 0.34 


®sc = 


0.00344 


P = 


0.14467 


Within populations 


186 11.14413 97.26 




0.02744 


P = 


0.00098 


(Madagascar/Maldives-Aceh/Java) 


Among groups 


1 -0.05540 -0.49 




-0.00490 


P = 


0.67742 


Among populations 


2 0.21697 1.92 


®sc = 


0.01910 


P = 


0.00293 


Within populations 


186 11.14413 98.57 


0 ST = 


0.01429 


P = 


0.00098 


(Madagascar/ Aceh-Maldives/ Java) 


Among groups 


1 -0.08441 -0.75 


0 cr = 


-0.00747 


P = 


1.00000 


Among populations 


2 0.23666 2.09 


®sc = 


0.02079 


P = 


0.00098 


Within populations 


186 11.14413 98.65 


0 ST = 


0.01348 


P = 


0.00098 



Central-eastern region (Maldives/ Aceh/Java) samples were pooled, suggested by population tree, 
haplotype network and ocean division. In addition, our preliminary analyses of pairwise @st comparisons 
of each data set found no heterogeneity as observed within samples from the central-eastern 
region (Maldives/Aceh/Java) (data not shown), although it was a limited Oct (Table 2, grouping 
Madagascar-Maldives/Aceh/Java). This supports the notion of a genetically homogeneous central-eastern 
region population that includes Maldives, Aceh and Java Sea, justified by combining the samples 
within the locality. 

Populations of P. penicillatus in the Indian Ocean region do not represent a genetically homogenous 
assemblage. Pairwise @st comparisons detected genetic structure between the central-eastern region 
(Maldives, Aceh, and Java) and other populations (Madagascar and Red Sea) (Table 3). A relatively 
small structure was found between the southwestern population (Madagascar) and central-eastern 
populations (Maldives/Aceh/Java) of P. penicillatus (@st value = 0.02629, p < 0.005; see Table 3). 
This pattern is in agreement with was previously observed from this region in reef fish 
Myripritis berndti [20] and amphidromous prawn Macrobrachium lar [21]. The prevailing ocean 
currents in the southwest Indian Ocean splits the coast of Madagascar into the southeast Madagascar 
current and the northeast Madagascar current, which may have limited genetic exchange with other 
regions [21,22]. 



Int. J. Mol. Sci. 2014, 15 



9248 



Table 3. Pairwise &st values (above the diagonal) and pairwise 0st ^-values (below the 
diagonal) of mitochondrial DNA control region among populations of Panulirus penicillatus. 



Population 


Northwestern 


Southwestern 


Central-Eastern Region 


Region (Red Sea) 


Region (Madagascar) 


(Maldives-Aceh-Java) 


Northwestern Region (Red Sea) 




0.51058 


0.53463 


Southwestern Region 
(Madagascar) 


0.00000 * 




0.02629 


Central-Eastern Region 
(Maldives- Aceh- Java) 


0.00000 * 


0.00098 * 





* Significant /?-value after Bonferroni correction (a = 0.05). 



The greatest significance is found in northwestern region (Red Sea) with supported by significant 
and high pairwise &st values (Table 3), UPGMA tree based on @ S t values showed a prominent genetic 
break between Red Sea and other populations (Figure 2) this result also corroborated with the 
distinction of the Red Sea cluster in the haplotype network (Figure 1), other studies on invertebrates, 
mud crab Scylla serrata [23,24], have shown a genetic isolation of Red Sea populations from the 
Indian Ocean. Oceanographic currents are likely to create barriers and impact routes and directions of 
larval dispersal [25-27]. More importantly, the lack of shared mtDNA haplotypes and large &st values 
suggest that there has been prolonged historical isolation between the Red Sea and others populations. 
In addition, high degree of endemism and biodiversity reef fishes were reported previously from 
the Red Sea, despite a geological history characterized by intermittent isolation and multiple 
salinity crises, mtDNA indicate that some widespread reef fish species originate in the Red Sea with 
subsequent contributing species to adjacent regions [28]. 

Further analysis on recruitment of P. penicillatus phyllosoma larvae may be necessary since genetic 
distinction in Red Sea populations may indicate the potential self-recruitment of phyllosoma larvae in 
this region. Although phyllosoma larvae are weak swimmers, they have the ability to maintain their 
vertical position in the water column [29]. Larval dispersal for spiny lobsters has been previously 
modeled. Griffin et al. [30], discovered that P. cygnus was a self-recruiting species, restricted to 
Western Australia. A more recent simulation study of P. argus found the average dispersal of P. argus 
in the Caribbean Sea to be possibly only 200-400 km, suggesting that larval behavior such as 
vertical migration, together with a retentive oceanographic environment, may increase the potential 
for self-recruitment [31]. 

No population structure was found in the localites of P. penicillatus inhabiting the central-eastern 
region of Indian Ocean: Nilandhoo atoll, Maldives; Aceh and Java Sea, Indonesia. The early life 
history of spiny lobsters consists of a drifting larval period adapted for a relatively long-term stay in 
the open ocean, extending from several months to more than a year, with many possibilities for 
dispersal through ocean currents [4,5]. Long larval life (8.3-9.4 months) of the P. penicillatus 
phyllosoma [32], coupled with the considerable impact of oceanic current systems or gyres (e.g., 
Equatorial currents and South Java current) (Figure 1) may explain the observed patterns. Furthermore, 
transportation of larvae likely occurred indirectly via stepping stones by the currents, and a small 
amount of larval dispersal over years might result in genetic homogeneity. Indeed, pelagic larval 
duration was considerably correlated with dispersal distance, but there were many exceptions; larval 
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behavior can play a crucial role in determining dispersal distance [33], while type of eggs and 
life-history parameters are also important predictors of connectivity in fishes [27,34]. 

Figure 2. UPGMA tree based on 0 ST values, illustrating the most probable geographical 
structure in the analysis of molecular variance (AMOVA). 

Red Sea, Jeddah, Saudi Arabia 

— I — Fort Dauphin, Madagascar 

r Nilandhoo atoll, Maldives 
Aceh, Indonesia 
Java Sea, Indonesia 

i 1 

0.05 

The neutrality of mutations in the mtDNA control region was rejected on the basis of Tajima's D 
and Fu' F tests. The significant negative values obtained in these P. penicillatus are often associated 
which experienced a population expansion during their geographic history [35]. This pattern has also 
been observed in other spiny lobster [12,36,37]. 

Although mtDNA has been proven to be a very useful marker, its use is not without complications. 
Due to the exclusively maternal inheritance mtDNA represents a marker strictly for historical process 
in females; they can potentially introgress between species [38,39], and symbiont-driven changes 
in mtDNA variation over space specifically in arthropods can occur [40]. Therefore, further studies 
incorporating supporting evidence from nuclear markers is required to obtain a more robust indication. 

3. Experimental Section 

3.1. Lobster Samples 

Adults of Panulirus penicillatus were collected from five localities in the Indian Ocean region 
including Red Sea, Jeddah, Saudi Arabia; Fort Dauphin, Madagascar; Nilandhoo atoll, Maldives; Aceh 
and Java Sea, Indonesia (Figure 1; Table 1). Lobsters were purchased from local commercial fishers 
and fully complied with local fisheries management and marine protected area controls. Tissues 
samples from walking legs (pereiopods) or abdomen were dissected on site, immediately fixed in 
70%-99% ethanol, and transferred to the laboratory. 

3.2. DNA Analysis 

About 50 mg of finely minced tissue sample was added to a 1.5 mL plastic test tube containing 
0.5 mL TNES-8 M urea buffer [41]. After adding 10 uL proteinase K, it was incubated at 38 °C for 
3 h and then genomic DNA was extracted using phenol-chloroform procedure and precipitated with 
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absolute ethanol according method modified in Imai et al. [42]. The mtDNA control region was 
amplified using the polymerase chain reaction (PCR). PCR primers previously designed for 
mitochondrial control region (Panulirusl2s: 5'-TATAGCAAGAATCAAACTATAG; and Penicillatus-R: 
5'-CATAGG(T/C)GTG(T/C)GAGGGAACAAAGTC) were used for amplification of the control 
region [3]. PCR amplifications were performed in 50 |xL reaction mixtures containing 1 jxL template 
DNA, 12.5 pmol of each primer and 25 |xL EmeraldAmp PCR Master Mix (2x Premix) (Takara Bio Inc., 
Shiga, Japan); the final volume of the reaction mixture was adjusted to 50 |xL with sterile water. 
Reactions were performed in a thermal cycler (Perkin Elmer GeneAmp PCR System 9700) under 
conditions of an initial denaturation step at 94 °C for 2 min, followed by 30 cycles of 30 s at 95 °C, 
30 s at 50 °C, and 1 min at 72 °C, with a final 7 min extension at 72 °C. PCR products were purified 
using a PCR Product Pre-sequencing kit (Exosap; USB Co., Cleveland, OH, USA). The cleaned 
products were sequenced on an ABI 3730x1 Genetic Analyzer (Applied Biosystems, Foster City, CA, 
USA) using a Big Dye Terminator Cycle Sequencing kit (ver. 3.1; Applied Biosystems). 

3.3. Genetic Data and Analyses 

Sequence data were aligned using ClustalX [43], with default alignment parameters and were 
checked manually for misalignments. The nucleotide compositions and numbers of variable sites were 
assessed with MEGA6 [44]. Haplotype and nucleotide diversity for each location were estimated using 
Arlequin (ver. 3.5) software [45]. AMOVA were performed to test the geographic divisions among 
population experimented with various groupings. This approach is a hierarchical approach that 
computes the proportion of variations among groups {Oct), the proportion of variation among 
population within groups ((£&), and the proportion of variation within populations (@st)- We ran 
preliminary analysis of pairwise 0st comparisons on each data sheet. The statistical significance of &st 
value was tested by 1000 permutations in Arlequin (ver. 3.5) software [45]. The Bonferroni test [46] was 
used to correct for multiple tests of the hypothesis that pairwise 0st statistics did not differ from zero. 
An UPGMA tree based on cp^r values was constructed using Neighbor in Phylip (ver. 3.6) [47,48]. The 
median-joining network [49], for the haplotypes was estimated using Network (ver. 4.611) [50]. Past 
demographic patterns were inferred by Tajima's D [51], and Fu's F [52], estimated using Arlequin 
(ver. 3.5). Both tests are commonly used to test neutrality; however, they can also be used to examine 
population growth because population expansion may result in rejection of the null hypothesis of 
neutrality (significant negative value). 

4. Conclusions 

The present study detected genetic structure in P. penicillatus populations located at the 
northwestern and southwestern edge of the Indian Ocean region with other populations. Relatively 
small structure was found between southwestern population and central-eastern, while vast majority 
structures were observed between the northwestern edge population and the other populations. Despite 
a lengthy larval period, oceanographic currents are likely to create barriers and impact routes and 
directions of larval dispersal in the northwestern edge population that may have limited genetic 
exchange with other regions. These results have implications for fisheries management in the region; 
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this should essentially provide information for ascertaining stock boundaries to better evaluate the 
long-term needs of spiny lobster management for this valuable species. 

At the scale of the central-eastern Indian Ocean, results suggested that P. penicillatus was belonging 
to a single unique panmictic population or at least several breeding grounds with significant exchange 
of genetic material. Long larval life (8.3-9.4 months) of the P. penicillatus coupled with the 
considerable impact of oceanic current systems or gyres (e.g., Equatorial currents and South Java 
current) may explain the observed patterns. Furthermore, transportation of larvae likely occurred 
indirectly via stepping stones by the currents, and a small amount of larval dispersal over some years 
might result in genetic homogeneity. 

Acknowledgments 

We thank the Merci Co., Ltd. (Shizuoka, Japan) for the help in sampling of spiny lobster from 
Madagascar. This study was supported in part by grants from the University of the Ryukyus grant for 
"International Research Hub Project for Climate Change and Coral Reef/Island Dynamics" and 
"the 2012 Nissui Research Fund" from Nippon Suisan Kaisha., Ltd. (Tokyo, Japan). 

Author Contributions 

Conceived and designed the experiments: M.F.A., H.I. Performed the experiments: M.F.A., H.I. 
Analyzed the data: M.F.A., H.I. Contributed to sampling materials/reagents: M.F.A., A., M.M., A.J.S., 

H. I. Wrote the paper: M.F.A., H.I. 

Conflicts of Interest 

The authors declare no conflict of interest. 
References 

I. Holthuis, L.B. Marine Lobsters of the World: An Annotated and Illustrated Catalogue of Species 
of Interest to Fisheries Known to Date; Food and Agriculture Organization of the United Nations: 
Rome, Italy, 1991; Volume 13. 

2. Pitcher, C.R. Spiny lobster. In Nearshore Marine Resources of the South Pacific. Information for 
Fisheries Development and Management; Wright, A., Hill, L., Eds.; IPS: Suva, Fiji/FFA: 
Honiara, Solomon/ICOD: Halifax, NS, Canada, 1993; pp. 539-607. 

3. Abdullah, M.F.; Chow, S.; Sakai, M.; Cheng, J.; Imai, H. Genetic diversity and population 
structure of pronghorn spiny lobster Panulirus penicillatus in the Pacific region. Pac. Sci. 2014, 
68, 197-211. 

4. Booth, J.D.; Phillips, B.F. Early life history of spiny lobster. Crustaceana 1994, 66, 271-294. 

5. Tolley, K.A.; Groeneveld, J.C.; Gopal, K.; Matthee, CA. Mitochondrial DNA panmixia in spiny 
lobster Palinurus gilchristi suggests a population expansion. Mar. Ecol. Prog. Ser. 2005, 297, 
225-231. 

6. Shaklee, J.B.; Bentzen, P. Genetic identification of stocks of marine fish and shellfish. 
Bull. Mar. Sci. 1998, 62, 589-621. 



Int. J. Mol. Sci. 2014, 15 



9252 



7. Diniz, F.M.; Maclean, N.; Ogawa, M.; Cintra, I.H.A.; Bentzen, P. The hypervariable domain of 
the mitochondrial control region in Atlantic spiny lobsters and its potential as a marker for 
investigating phylogeographic structuring. Mar. Biotechnol. 2005, 7, 462-473. 

8. Avise, J.C. Molecular Markers: Natural History and Evolution; Chapman and Hall: London, 
UK, 1994. 

9. Inoue, N.; Watanabe, H.; Kojima, S.; Sekiguchi, H. Population structure of Japanese spiny lobster 
Panulirus japonicus inferred by nucleotide sequence analysis of mitochondrial COI gene. 
Fish. Sci. 2007, 73, 550-556. 

10. Garcia-Rodriguez, F.J.; Perez-Enriquez, R. Lack of genetic differentiation of blue spiny lobster 
Panulirus inflatus along the Pacific coast of Mexico inferred from the mtDNA sequences. 
Mar. Ecol. Prog. Ser. 2008, 361, 203-212. 

11. Sarver, S.; Silberman, J.D.; Walsh, P.J. Mitochondrial DNA sequence evidence supporting 
the recognition of two subspecies or species of the Florida spiny lobster Panulirus argus. 
J. Crust. Biol. 1998, 18, 177-186. 

12. Farhadi, A.; Farhamand, H.; Nematollahi, M.A.; Jeffs, A.; Lavery, S.D. Mitochondrial DNA 
population structure of the scalloped lobster Panulirus homarus (Linnaeus 1758) from the West 
Indian Ocean. ICES J. Mar. Sci. 2013, 70, 1491-1498. 

13. Chow, S.; Jeff, A.; Miyaeke, Y.; Konishi, K.; Okazaki, M.; Suzuki, N.; Kimura, S.; Abdullah, M.F.; 
Imai, FL; Sasaki, M.; et al. Genetic isolation between the western and eastern Pacific populations 
of pronghorn spiny lobster Panulirus penicillatus. PLoS One 2011, 6, e29280. 

14. Schott, F.; McCreary, J.P. The monsoon circulation of the Indian Ocean. Prog. Oceanogr. 2001, 
51, 1-123. 

15. McMillen- Jackson, A.L.; Bert, T.M. Genetic diversity in the mtDNA control region and 
population structure in the pink shrimp Farfantepenaeus duorarum. J. Crust. Biol. 2004, 24, 
101-109. 

16. McMillen-Jackson, A.L.; Bert, T.M. Disparate patterns of population genetic structure and 
population history in two sympatric penaeid shrimp species (Farfantepenaeus aztecus and 
Litopenaeus setiferus) in the eastern United States. Mol. Ecol. 2003, 12, 2895-2905. 

17. Imai, H.; Hanamura, Y.; Cheng, J.H. Genetic and morphological differentiation in the Sakura 
shrimp (Sergia lucens) between Japanese and Taiwanese populations. Contrib. Zool. 2013, 82, 
123-130. 

18. Heyer, E.; Zietkiewicz, E.; Rochowski, A.; Yotova, V.; Puymirat, J.; Labuda, D. Phylogenetic and 
familial estimates of mitochondrial substitution rates: study of control region mutations in 
deep-rooting pedigrees. Am. J. Hum. Genet. 2001, 69, 1 113-1 126. 

19. Bird, C.E.; Smouse, P.E.; Karl, S.A.; Toonen, R.J. Detecting and measuring genetic 
differentiation. In Crustacean Issues: Phylogeography and Population Genetics in Crustacea; 
Koenemann, S., Held, C, Schubart, C, Eds.; CRC Press: Boca Raton, FL, USA, 2011; pp. 31-55. 

20. Muths, D.; Tessier, E.; Gouws, G.; Craig, M.; Mwale, M.; Mwaluma, J.; Mwandya, A.; Bourjea, J. 
Restricted dispersal of the reef fish Myripristis berndti at the scale of the SW Indian Ocean. 
Mar. Ecol. Prog. Ser. 2011, 443, 167-180. 



Int. J. Mol. Sci. 2014, 15 



9253 



21. Castelin, M.; Feutry, P.; Hautecoeur, M.; Marquet, G.; Wowor, D.; Zimmeraiann, G.; Keith, P. New 
insight on population genetic connectivity of widespread amphidromous prawn Macrobrachium lar 
(Fabricius, 1798) (Crustacea: Decapoda: Palaemonidae). Mar. Biol. 2013, 160, 1395-1406. 

22. Chapman, P.; di Marco, S.F.; Davis, R.E.; Coward, A.C. Flow at intermediate depths around 
Madagascar based on ALACE float trajectories. Deep Sea Res. Part II Top. Stud. Oceanogr. 
2003, 50, 1957-1986. 

23. Gopurenko, D.; Hughes, J.M.; Keenan, CP. Mitochondrial DNA evidence for rapid colonisation 
of the Indo-West Pacific by the mudcrab Scylla serrata. Mar. Biol. 1999, 134, 227-233. 

24. Fratini, S.; Vannini, M. Genetic differentiation in the mud crab Scylla serrata (Decapoda: 
Portunidae) within the Indian Ocean. J. Exp. Mar. Biol. Ecol. 2002, 272, 103-1 16. 

25. Baums, I.B.; Paris, C.B.; Cherubin, L.M. A biooceanographic filter to larval dispersal in a 
reef-building coral. Limnol. Oceanogr. 2006, 51, 1969-1981. 

26. Treml, E.A.; Halpin, P.N.; Urban, D.L.; Pratson, L.F. Modeling population connectivity by ocean 
currents, a graph-theoretic approach for marine conservation. Landsc. Ecol. 2008, 23, 19-36. 

27. Riginos, C; Douglas, K.E.; Jin, Y.; Shanahan, D.F.; Treml, EA. Effects of geography and life 
history traits on genetic differentiation in benthic marine fishes. Ecography 2011, 34, 566-575. 

28. DiBattista, J.D.; Berumen, M.L.; Gaither, M.R.; Rocha, LA.; Eble, J.A.; Choat, J.H.; Craig, M.T.; 
Skillings, D.J.; Bowen, B.W. After continents divide: comparative phylogeography of reef fishes 
from the Red Sea and Indian Ocean. J. Biogeogr. 2013, 40, 1 170-1 181. 

29. Johnson, M.W. On the dispersal of lobster larvae into the East Pacific Barrier (Decapoda, 
Palinuridea). Fish. Bull. 1974, 72, 639-647. 

30. Griffin, DA.; Wilkin, J.L.; Chubb, C.F.; Pearce, A.F.; Caputi, N. Ocean currents and the larval phase 
of Australian western rock lobster, Panulirus cygnus. Mar. Freshw. Res. 2001, 52, 1 187-1 199. 

31. Butler, M.J., IV; Paris, C.B.; Goldstein, J.S.; Matsuda, H.; Cowen, R.K. Behavior constrains the 
dispersal of long-lived spiny lobster larvae. Mar. Ecol. Prog. Ser. 2011, 422, 223-237. 

32. Matsuda, H.; Takenouchi, T.; Goldstein, J.S. The complete larval development of the pronghorn 
spiny lobster Panulirus penicillatus (Decapoda: Palinuridae) in culture. J. Crust. Biol. 2006, 26, 
579-600. 

33. Shanks, A.L. Pelagic larval duration and dispersal distance revisited. Biol. Bull. 2009, 216, 
373-385. 

34. Bradbury, I.R.; Laurel, B.; Snelgrove, P.V.R.; Bentzen, P.; Campana, S.E. Global patterns in 
marine dispersal estimates: The influence of geography, taxonomic category and life history. 
Proc. R. Soc. B 2008, 275, 1803-1809. 

35. Rogers, A.R.; Harpending, H. Population growth waves in the distribution of pairwise genetic 
differences. Mol. Biol. Evol. 1992, 9, 552-569. 

36. Palero, F.; Abello, P.; Macpherson, E.; Gristina, M.; Pascual, M. Phylogeography of the European 
spiny lobster (Palinurus elephas): Influence of current oceanographical features and historical 
processes. Mol. Phylogenet. Evol. 2008, 48, 708-717. 

37. Babbucci, M.; Buccoli, S.; Cau, A.; Cannas, R.; Goni, R.; Diaz, D.; Marcato, S.; Zane, L.; 
Patarnello, T. Population structure, demographic history, and selective processes: 
Contrasting evidences from mitochondrial and nuclear markers in the European spiny lobster 
Palinurus elephas (Fabricius, 1787). Mol. Phylogenet. Evol. 2010, 56, 1040-1050. 



Int. J. Mol. Sci. 2014, 15 



9254 



38. Ballard, J.W.; Whitlock, M.C. The incomplete natural history of mitochondria. Mol. Ecol. 2004, 
13, 729-744. 

39. Jiggins, F.M. Male-killing Wolbachia and mitochondrial DNA: Selective sweeps, hybrid 
introgression and parasite population dynamics. Genetics 2003, 164, 5-12. 

40. Hurst, G.D.D.; Jiggins, F.M. Problems with mitochondrial DNA as a marker in population, 
phylogeographic and phylogenetic studies: The effect of inherited symbionts. Proc. R. Soc. B 
2005, 272, 1525-1534. 

41. Asahida, T.; Kobayashi, T.; Saitoh, K.; Nakayama, I. Tissue preservation and total DNA 
extraction from fish stored at ambient temperature using buffers containing high concentration of 
urea. Fish. Sci. 1996, 62, 727-730. 

42. Imai, H.; Cheng, J.-H.; Hamasaki, K.; Numachi, K.-I. Identification of four mud crab species 
(genus Scylla) using ITS-1 and 16S rDNA markers. Aquat. Living Resour. 2004, 17, 31-34. 

43. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The ClustalX 
windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis 
tools. Nucleic Acids. Res. 1997, 25, 4876-4882. 

44. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary 
Genetics Analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725-2729. 

45. Excoffier, L.; Laval, G.; Schneider, S. Arlequin ver. 3.0: An integrated software package for 
population genetics data analysis. Evol. Bioinform. Online 2005, 1, 47-50. 

46. Rice, W.R. Analyzing tables of statistical tests. Evolution 1989, 43, 223-225. 

47. Felsenstein, J. PHYLIP-Phylogeny Inference Package. Cladistics 1989, 5, 164-166. 

48. Felsenstein, J. PHYLIP Home Page. Available online: http://evolution.genetics.washington.edu/ 
phylip.html (accessed on 21 May 2014). 

49. Bandelt, H.J.; Forster, P.; Rohl, A. Median-joining networks for inferring intraspecific 
phylogenies. Mol. Biol. Evol. 1999, 16, 37-48. 

50. Fluxus-engineering.com. Available online: http://www.fluxus-engineering.com (accessed on 21 
February 2014). 

51. Tajima, F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. 
Genetics 1989, 123, 585-595. 

52. Fu, Y.-X. Statistical tests of neutrality of mutations against population growth, hitchhiking and 
backgroud selection. Genetics 1997, 147, 915-925. 



© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http://creativecommons.Org/licenses/by/3.0/). 



